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Abstract: We study the feasibility of nonlocally compensating for polar- 
ization mode dispersion (PMD), when polarization entangled photons are 
distributed in fiber-optic channels. We quantify the effectiveness of nonlocal 
compensation while taking into account the possibility that entanglement is 
generated through the use of a pulsed optical pump signal. 
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1. Introduction 



Applications of quantum information protocols often rely on the ability of distributing quan- 
tum entanglement among distant users, a task for which the installed fiber-optic infrastructure 
is the most natural candidate. Two types of photon-entanglement are particularly popular; the 
entanglement of polarization 11], and time-bin entanglement 13. While the latter scheme has an 
important advantage of being more tolerant to fiber propagation effects, the former is attractive 
due to the ease with which photon polarizations can be manipulated by using standard opti- 
cal instrumentation. The greatest disadvantage of the polarization-entanglement scheme stems 
from its inherent sensitivity to the presence of polarization decoherence mechanisms O, the 
most important of which is polarization mode dispersion (PMD) 0. This phenomenon, which 
is generated by the random residual birefringence of optical fibers, introduces a frequency de- 
pendent polarization rotation of the propagating signals, a process whose effect on polarization- 
entanglement could potentially be devastating |5 |. As the consequences of PMD can be very 
severe also in classical high-data-rate optical communications systems, many schemes for its 
compensation have been proposed IS). The same compensating schemes could in principle 
be applied independently to each fiber in an entanglement distribution system. On the other 
hand, quantum entanglement has an intriguing non-classical property that allows approaching 
the compensation problem in a somewhat unusual form. The combined effect of PMD in each 
fiber connecting the two receiving stations (conventionally referred to as Alice and Bob) can in 
principle be eliminated by introducing a compensator in only one of the two optical paths, for 
example in the path of the photon destined for Bob. This concept, to which we refer as nonlo- 
cal PMD compensation, is similar to Franson's nonlocal compensation of chromatic dispersion 
Qm, and it is closely related to the concept of decoherence-free subspaces, demonstrated by 
Kwiat 121 and Altepeter [ 10| in optically birefringent crystals. 

In this paper we consider the feasibility of nonlocal PMD compensation in fiber-based entan- 
glement distribution systems. In particular, we concentrate on the case in which the entangled 
photons are generated in a nonlinear medium pumped by a pulsed optical source flTI . The 
regime of pulsed pumping is of particular importance to applications as many schemes for en- 
tanglement generation in optical fibers use pumps that consist of optical pulses lfT2l [131 [T4l . 
In what follows, we explore the relation between the pump-pulse duration (or bandwidth), the 
bandwidth of the entangled photons and the quality of compensation. 

2. PMD as a polarization decoherence mechanism 

In its simplest manifestation, valid in the case of sufficiently narrow-band signals, the effect of 
PMD is to cause an incident optical pulse to be split into two orthogonally polarized compo- 
nents that are delayed relative to each other The polarization states of these two components 
are known as the principal states of polarization (PSP) and the delay between them is called the 
differential group delay (DGD) Q. This description is commonly referred to as the first-order 
picture, and within its range of validity, the entire phenomenon of PMD is similar to that of pure 
birefringence. An important difference is that in the case of PMD both the DGD and the PSPs 
vary randomly with frequency and time 11511 . The first order picture is valid as long as the band- 
width characterizing the spectral evolution of the PMD parameters remains much larger than 
the bandwidth of the signals transmitted through the fiber. The former is given approximately 
by the inverse of the average DGD value ifTSl . In order to preserve the clarity of the arguments 
and results presented in this paper, we conduct the discussions that follow for the case in which 
the first-order picture applies. Nonetheless, the principles underpinning our discussions can be 
readily (albeit in an very cumbersome form) expanded to the most general case. 

To illustrate the way in which PMD causes decoherence and then to introduce the concept 
of nonlocal compensation, we start by discussing the case of a perfectly stationary stream of 
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Fig. 1. Nonlocal PMD compensation set-up. The compensator in Bob's path contains two 
stages. The first compensates Bob's PMD, the second compensates for Alice's PMD. 

entangled photons, generated by a continuous wave (CW) pump. Assume first that PMD is 
present only in Alice's arm. In this case her photon can either be detected earlier, or later than 
the photon registered by Bob's detector, thereby disclosing that its state of polarization has col- 
lapsed to either the fast, or the slow PSP, respectively. Since the receivers are insensitive to the 
photons' arrival times, this information is lost to the environment, resulting in decoherence and 
in a reduction in the degree of polarization entanglement. We may now introduce the concept 
of nonlocal compensation. Assume that Bob includes in his optical path a PMD element whose 
properties over the signal spectrum are identical to those of Alice's. His photon-being polar- 
ization entangled with Alice's-will necessarily collapse into the same PSP (slow, or fast) as 
Alice's photon, so that the two photons will be detected simultaneously. In this case, the times 
of arrival disclose no information on whether the photons passed through the slow, or the fast 
PSP and decoherence is avoided. The situation changes drastically when a pulsed pump is used 
in order to generate the entangled photons. In this case, the pump imposes an independent time 
reference, which enables one to assess whether the photons are aligned with the slow, or the fast 
PSP, based on their time of arrival relative to the pump's time reference. It is therefore appar- 
ent that in this scenario, decoherence is unavoidable and the effectiveness of the compensation 
process deteriorates. In the section that follows we provide a formal description of this regime 
and quantitatively assess the compensation efficiency. 

3. Theory and results 

Figure 1 shows a schematic description of the considered nonlocal compensation scheme. Bob's 
compensator is illustrated as a two-stage device, where the first stage is responsible for compen- 
sating Bob's PMD, whereas the second stage is the nonlocal compensator which compensates 
for the PMD present in Alice's fiber. While this description allows us to ignore the PMD that 
is present in Bob's arm in what follows, we point out that in a realistic implementation both 
stages of compensation can be realized in a single device. We consider a maximally polariza- 
tion entangled two-photon state, generated by parametric down conversion in a nonlinear 
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The functions Ha{co) and Hi,{co) are the transfer functions of two frequency filters positioned 
along Alice and Bob's optical paths, respectively, and Ep{co) is the spectrum of the pump signal. 
We have assumed in Eq. ([U that the phase matching condition is uniformly satisfied within the 
bandwidth of the two filters, so that the spectral dependence of phase matching is negligible. 
The terms u^, u^, and denote the Jones vectors describing the polarization states of Alice 
and Bob's photons, with the primes denoting orthogonality in polarization space, namely m^, • 

'The same expression with small modifications to the waveform part applies to the case where entanglement is 
generated in a nonlinear optical medium 

possible phase difference between the two polarization terms, which often follows from the experimental pro- 
cedure of photon generation |12| , is immaterial to our analysis and is therefore omitted. 




= Mi ■ = 0- The frequency dependent part in Eq. ^ can be conveniently reexpressed as 

^ydCOa^ d(0bH*{(0a)Hl{(0b)Ep{(0a + (0h)\(0a,(0h) = j dt,, j dti,g{ta,tb)\ta,tb) (2) 

where g{ta,tb) = / df /!*,(? — ti,)hl{t -'tb)Ep{t), with ha{t) and hb{t) being the impulse responses 
of the two filters and where Ep{t) represents the complex envelope of the pump signal. As 
is apparent from Eq. (|2]i, |g(f„,ff,)p is proportional to the probability density that Alice and 
Bob's photons arrive at the detectors at the times f„ and tb, respectively. Normalization of the 
overall quantum state requires that g{ta,tb) satisfies / dtadtb\g{ta,tb)\^ = 1. In the case of CW 
pumping, Ep is a constant and g{ta,tb) reduces to be a function of the time difference T = r„ — tb, 
thereby making the entangled photon stream stationary ||7][8l|9l. In this case the normalization 
condition must be modified to / dT|g(T) p = ' , with T being the photo-detector's integration 
window. In order to abbreviate the notation, we denote the polarization part of the state-vector 
by I w) = 2^'/^[|m^,,m^) + |m^i,mJ,)], and by \g{ta,tb)) the waveform related part given by Eq. 
(|2]i. This way, the representation of the overall state simplifies to | y/) = \g{ta,tb)) <E) \Wp)- 

We now consider the case in which Bob introduces a second PMD element into his optical 
path in order to compensate for PMD that is suffered by Alice's photon (stage #2 in Fig. 1). In 
order to conveniently describe the effect of PMD on the system, we note that the two-photon 
polarization state \\(fp) maintains its form in a subspace of representation bases. In particular, 
if one denotes the slow PSP in Alice's arm by i^,, one can always find a state of polarization 
Sj, such that \\j/p) is expressed in the {i;^,,^^} basis as \\j/p) = 2^^^^[\s^,Sb) + \s!cn^'b)]- ^'^^ 
apparent that in order to perform nonlocal PMD compensation, stage #2 of Bob's compensator 
must contain a PMD element, whose slow PSP is aligned with i^. Then the quantum state of 
the photons received by Alice and Bob can be conveniently expressed as 

Irout) - -^Isita - ra/2,tb - Tb/2)) ® \s^„S^) + + ^J2,tb + Tb/2)) ® 1*^,^;), (3) 

where Ta is Alice's DGD and Tb is the DGD in stage #2 of Bob's compensator It is ev- 
ident in Eq. (O how the photons' states of polarization become coupled to their times of 
arrival. The degree of entanglement can be assessed by calculating the concurrence ITtI 
from the density matrix that corresponds to jvAjut) after tracing it over the time modes, i.e. 
p = J dtadtb{ta,tb\Wout) {Wom\ta,tb)- Following some algebra, the concurrence can be presented 
in the following form 
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where normalization of g{ta,tb) ensures that in the absence of PMD (when T„ = T/, = 0), the 
concurrence is equal to unity. The existence of a decoherence-free subspace demonstrated in 
ll9l. |[T0l follows from reducing the bandwidth of the pump signal to zero and by setting the 
DGD values T„ and Tb to be equal. In this case, as we mentioned earlier, normalization of g{T) 
implies that C = 1 and the degree of entanglement is fully preserved. Note that this is a direct 
consequence of the fact that both photons are either delayed (if they are in the \s^i,s^) state), 
or advanced (if they are in the \s^,,s^b) state). As for a continuous wave pump there is no abso- 
lute time reference that would indicate whether they are late, or early, no information on their 
polarization states can be extracted from their times of arrival, thereby implying the absence 
of decoherence. Conversely, when the the pump signal consists of pulses of finite duration, the 
photons' times of arrival relative to the pump disclose information on their polarizations, even 
if they reach the detectors simultaneously. In this scenario decoherence cannot be avoided. 



We now quantify the above principles by considering an example consisting of Gaussian 
shaped frequency filters and of a pump signal whose complex envelope is a Gaussian pulse. 
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, where we use 



Namely, Ep{o)) °^ exp (-(0^/4fip) and Ha.b{(o) °= exp 

to denote the root-mean-square bandwidth of the optical pump |£'((b)|^ and the term Z?„ f, de- 
notes the root-mean-square bandwidths of Alice's, or Bob's filters |//„^(a))p. The frequency 
axis was chosen such that the pump frequency is equal to and the central frequencies of Alice 
and Bob's filters are symmetrically offset from it by iAH. In this condition Eq. (|4| reduces to 
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The first term on the right-hand-side of (|5]i accounts for the loss of entanglement due to a 
difference in the DGD values, whereas the second term represents the main effect of the pump 
bandwidth. Equation (|5]l reveals the relevant regimes of operation. When the pump bandwidth 
is small relative to the bandwidths Ba and B^, of Alice and Bob's photons, the second term in 
Eq. (|5]l reduces to unity. Then the overall concurrence can be made equal to 1 if the DGD in 
Bob's compensator is equal to AUce's DGD, i.e if ~ To. In the opposite Hmit, when Bp is 
greater than Ba and Bj,, the concurrence assumes a value that depends on the magnitude of the 
DGD. Provided that the DGD is not significantly larger than the inverse bandwidth of Alice 
and Bob's filters, the concurrence in this case reduces to C = exp [—[B^tI+b'^'zI) /2] , a value 
that is independent of Bp. The highest concurrence that can be achieved by nonlocal PMD 
compensation is given by 
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and it is obtained when the compensators DGD is given by T°P' = T„ / (1 +Bl/Bl),a value that 
naturally reduces to T^'" = T„ when the pump bandwidth is set to Bp = 0. On the other hand, 
T^P' approaches 0, when Bp ^ Z?„. Since in most situations of practical interest. Bp < B„, effec- 
tive compensation takes place as long as the DGD value to be compensated is small relative to 
Bp i.e. relative to the pump-pulse duration. Notice also that for Gaussian waveforms neither 
the optimal value of T/, needed for compensation, nor the compensated concurrence depend on 
the bandwidth of Bob's optical filter. Finally, we note that the sensitivity to small errors in the 
setting of the compensator's DGD is very weakly dependent upon the pump bandwidth, and 
thus for all practical purposes, the concurrence can be approximated as C ~ C°P' ( 1 — 5t^/Tq) , 



with Xb = If + 5t and where Tq ^ 



IBp^ 



Interestingly, in a short (^lOOKm) link, of the kind 



used in repeaterless quantum communications, even if the instantaneous DGD reaches values 
of the order of few picoseconds, compensation allows maintaining the concurrence above 0.95 
for typical pump and channel bandwidths of lOOGHz {Bp = B„ = 27rx lOOGHz). For reference, 
note that acceptable quantum bit error rates in typical applications of quantum protocols usually 
translate into concurrence values in excess of 0.6. The above analysis assesses the maximum 
achievable concurrence as a function of the pump signal's bandwidth, and provides the DGD 
value that needs to be assigned to the compensator for optimal performance. While these results 
are analytical in the case of Gaussian waveforms and filters, they provide a fairly good estimate 
for the situation in the more realistic scenarios. To see that, we numerically evaluate the con- 
currence in a scenario that is consistent with the main features of a commercially available 
entanglement generation system lfT2l[T3]| . which uses typical telecom filters for both the pump 
and the signals. As we found in measurements |I3|, the spectra of these filters can be accurately 
described by a third-order super-Gaussian profile. In all the numerical results presented in what 
follows, we assume that Alice's and Bob's filters have identical bandwidths Ba = Bi, = B. 



Fig. 2. (a) The concurrence as a function of tire tire DGD of Bob's compensator normalized 
by the DGD in Alice's optical path. The red and blue curves correspond to pump band- 
widths of Bp = 0. 1 ' and Bp = T^^, respectively, (b) The best achievable concurrence af- 
ter PMD compensation, as a function of the pump bandwidth Bp normalized by ' , where 
the green and black curves correspond to photon bandwidths of 6 = 2t^^ and B = z^^, 
respectively. The solid curves are drawn for typical telecom filter-shapes (super-Gaussian), 
whereas the dashed curves show the analytical expressions in the Gaussian case. 

In Fig. 2a we plot the concurrence versus the DGD of Bob's compensator, for the case where 
TaB = 1 . The results corresponding to the super-Gaussian waveforms are represented by the 
solid curves, whereas the dashed curves illustrate Eq. (|5]l, for comparison. The two sets of colors 
correspond to two values of the pump bandwidth; the case of Bp — 0.1/ Ta is shown in blue, and 
Bp = 1/ T„ is shown in red. In the case of Bp = 0.1 /t„, namely when the DGD is much smaller 
than the time duration of the pump pulses, perfect compensation can be achieved when ~ Ta- 
in the other case, when Bp = 1 / T„, so that the pump duration is comparable to the compensated 
DGD, the curve peaks at a smaller value of and the optimal concurrence is visibly lower than 
unity. While the detailed features of the dashed and solid curves differ somewhat, the general 
properties are similar, as expected. Figure 2b shows the maximum achievable concurrence as 
a function of the pump bandwidth Bp, represented in units of T^^. In this case, the two sets 
of curves correspond to two different bandwidths of Alice and Bob's filters; the blue curves 
correspond to B — 1/To, whereas the red curves correspond to B = 2/t„. The width of the 
compensated concurrence functions in Fig. 2b gives an idea on the range of pump bandwidths 
that allows effective compensation. As can be seen from Eq. (|6]l, high quality compensation 
occurs when Bp is smaller that the inverse DGD that is to be compensated, i.e. Bp < z^K 
Beyond that value, the remaining entanglement depends on the bandwidth of the entangled 
photons. The smaller the value of B, the higher the value to which the concurrence reduces 
when the bandwidth of the pump pulses is increased. This is not surprising as in the limiting 
case, where B <C the temporal uncertainty in aiTival times of the entangled photons is so 
long, that PMD has nearly no effect on it, and the degree of entanglement remains close to 
unity. In this limit, however, PMD compensation is of little relevance. 

4. Conclusions 

To conclude, we have studied the concept of nonlocal PMD compensation in the context of 
distributing pairs of polarization entangled photons over a fiber-optic link. While perfect com- 
pensation can be achieved when the stream of photons is stationary, the quality of compensation 
reduces rapidly when entanglement generation schemes that rely on pulsed pumping are used. 
For the special case of Gaussian waveforms, we analytically found the optimal compensator 
setting and the resulting quality of compensation. Numerical estimation of compensation in a 
system with realistic telecom filter shapes produced results close to the analytical ones. 



